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Abstract 

We compare the lower edge spectral fluctuations of the staggered 
lattice Dirac operator for the Schwinger model with the predictions 
of chiral Random Matrix Theory (chRMT). We verify their range of 
applicability, checking in particular the role of non-trivial topological 
sectors and the flavor symmetry of the staggered fermions for finite 
lattice spacing. Approaching the continuum limit we indeed find clear 
signals for topological modes in the eigenvalue spectrum. These find- 
ings indicate problems in the verification of the chRMT predictions. 
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1 Introduction 



Leutwyler and Smilga's sum rules pQ for the eigenvalues of the chiral Dirac 
operator for finite volume QCD suggest that the statistical fluctuations of the 
spectrum of this operator are universal in the infrared region; these should 
be described by a chiral Random Matrix Theory (chRMT) (For a recent 
review on Random Matrix Theory in general cf. [3 .) Universality has been 
proved extensively in chRMT; for the case of QCD, see [3]. Explicit lattice 
calculations may provide an answer to the question, whether QCD lies really 
inside one of the universality classes of chRMT. 

Chiral symmetry, which is the main ingredient here, may be effectively 
realized on the lattice by the so-called "overlap fermions" [5: and by the 
Fixed-Point (FP) fermion action jB] , both satisfying the Ginsparg- Wilson [7] 
condition. Ginsparg- Wilson fermions - in spite of their theoretical appeal - 
present however non-trivial technical difficulties in a dynamical set-up (for 
first experiences in this context, see [BHH!)- 

Kogut-Susskind "staggered" fermions, which are the object of this letter, 
realize chirality only partially, but are also much simpler, and still represent 
the typical framework for lattice QCD calculations when chirality is relevant. 
However, this is the main conclusion of this letter, comparison with chRMT is 
in this case not straightforward and subject to additional restrictions (besides 
the usual ones well-known in the literature) applying as long as the lattice 
cut-off is finite. 

A class of problems comes from topology: when comparing lattice eigen- 
value spectra with chRMT predictions, the validity of the index theorem [TO] 
for the Dirac operator is assumed, while the staggered Dirac operator has no 
exact zero modes and a sound fermionic definition of the topological charge 
is not possible. Another class of problems is due to the flavor structure of 
staggered fermions, which is different at finite cut-off and in the continuum 
limit. Each copy of staggered fermions describes in the continuum Nj = 2 d ^ 2 
degenerate massless physical fermions; at finite cut-off the full chiral sym- 
metry is broken into an abelian subgroup resembling the non-degenerate 
(Nf = 1) situation, with the relevant difference that the symmetry is here 
anomaly-free. 

We study these issues in the simplified framework of the Schwinger model 
in the quenched and the dynamical situation. A particularity of two dimen- 
sions is that a continuous symmetry cannot be broken by the vacuum [TT] . 
Since the chiral symmetry at finite cut-off is non-anomalous though abelian, 
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the usual fermion condensate vanishes even if staggered fermions act as just 
one flavor when comparing Dirac eigenvalue spectra with chRMT. With van- 
ishing condensate, a different kind of universality should apply in the infrared 
part of the spectrum [T2] . 

In a Monte Carlo simulation, all this applies of course only for the un- 
quenched (dynamical) data. For the "quenched" model (corresponding to 
Nf = 0) the usual chRMT universality is expected to hold |B]. In this case, 
however, different complications come into play when the thermodynamic 
limit of the theory is taken [T3*] . Some results concerning the quenched data 
have been discussed in [14j. 

2 Predictions from Random Matrix Theory 

The lower edge of the spectrum of the Dirac operator (i.e. its infra-red region) 
is expected to be universal at the scale of the level spacing AA (microscopic 
scale) ~ 1/V". The 'microscopic' scaling variable is z = 7 S A, where 



and p(A) denotes the associated spectral density per unit volume. When 
comparing the statistical properties of the Dirac operator with chRMT, E is 
an external scale parameter related to the dynamics of the physical system 
under study. The Banks-Casher formula jT^I implies S = —(ipip). 

Three classes of universality are predicted by chRMT [16J, corresponding 
to (chiral) orthogonal, unitary and symplectic ensembles (chOE, chUE and 
chSE respectively). The universal properties considered in this letter are the 
microscopic spectral density 



and the probability distribution of the smallest eigenvalue P(A m i n ). These 
distributions depend on the topological charge v of the background gauge 
configuration and on the number of flavors Nf of quarks included in the 
dynamics. We also checked the level spacing statistics [T3] . 



£ = 7T lim lim 



(1) 




(2) 
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3 Statistical properties of the spectrum 



Staggered fermions partially realize chiral symmetry on the lattice. One 
species of staggered fermion describes Nf = 2 d ^ 2 Dirac fermions in the contin- 
uum limit. For finite lattice spacing the chiral symmetry £7y(2 d / 2 ) x Ua^^ 2 ) 
is however broken down to a UyiX) x ^7a(1) non-anomalous sub-symmetry, 
which accounts for the reflection-invariance for the purely imaginary spec- 
trum. 

In chRMT, the number of dynamical fermions Nf and the topological 
charge v of the background configuration enter as external parameters. Our 
(compact) gauge configurations were sampled in the quenched set-up, i.e. 
according to the pure gauge measure (given by the standard plaquette ac- 
tion). This corresponds in chRMT to the case Nf = 0. Assuming no extra 
symmetry for the Dirac operator (besides the one mentioned above) in our 
two-dimensional set-up, the statistical fluctuations of the spectrum in the 
infrared region are expected ^Bj to be described by chUE; this was con- 
firmed by previous results [H] with the Fixed Point and Neuberger's overlap 
action, where however, deviations (i.e. a distribution resembling chSE) were 
observed whenever the physical volume was too small. 

In order to explore the dynamics of fermions, we then re-weighted the 
quenched configurations with the fermion determinant. For d = 4 this proce- 
dure usually introduces large statistical fluctuations; in our two-dimensional 
model these fluctuations did not appear to be too harmful |5]. Since the 
chiral symmetry of lattice staggered fermions is abelian, one expects ^7] 
that the dynamical situation is described by a chRMT with Nf = 1. The 
two-dimensional situation is however peculiar due to the Mermin-Wagner- 
Coleman theorem stating that a continuous symmetry cannot be broken by 
the vacuum; only an anomaly can do that. Since the operator ip(x)ip{x) is 
not invariant under the residual (non-anomalous) abelian chiral symmetry of 
staggered fermions, its expectation value vanishes [TH]. (One can argue [TH] 
that the simplest observable signaling the breaking of the anomalous axial 
symmetry is in this context (uudd).) With vanishing fermion condensate, 
the universality class of chUE does not apply anymore, and a new kind of 
universality holds |12j . 

The validity of the Nf = 1 chRMT predictions is restricted. In the contin- 
uum limit the irrelevant terms of the action breaking the full symmetry vanish 
and a transition to Nf = 2 is due to occur with doubling of the spectrum. In 
this limit, two degenerate fermions described by two-component spinors are 



4 



the effective degrees of freedom, the two-flavor (massless) Schwinger model 
being the underlying continuum theory [TK] . 

RMT predictions depend also on the topological charge of the background 
gauge configuration. In the continuum the topological charge is related to 
the index of the Dirac operator JUj (the number of zero modes counted 
according to their chirality). On the lattice, a consistent definition of the 
topological charge in terms of the index of the Dirac operator is possible 
only for Ginsparg- Wilson (i.e. Overlap and FP) fermions [20] • For staggered 
fermions, away from the continuum limit, the Dirac operator has no exact 
zero modes, and the naive expectation is that all configurations show spectral 
distributions described by chRMT predictions for the trivial sector [T7| I21j. 

However, for lattices that are fine enough one should attain the contin- 
uum situation. We will explore this hypothetical scenario and check possible 
systematic effects when comparing the spectrum of the Dirac operator with 
chRMT predictions. Lacking (in this case) a consistent fermionic definition 
of the topological charge v, we rely on the geometric definition of the gauge 
configurations. In the analysis we distinguish between the sectors with trivial 
[y = 0) and non-trivial topological charge. 

3.1 The simulation 

We considered statistically independent gauge configurations for lattices of 
size 16 2 (p = 2, 4, 6; 10000 configurations); 24 2 {(3 = 2, 4, 6, 8 ; 5000 
configurations) and 32 2 ( j3 — 2, 4, 6; 5000 configurations). Even-odd pre- 
conditioning and standard LAPACK routines were used to diagonalize the 
staggered Dirac operator. 

The trivial sector. Fig. ^shows the distribution of the smallest eigenvalue 
-P(Amin) for some lattice volumes. The histograms are obtained by choosing 
configurations with geometrical topological charge v = 0; this reduces the 
statistics to 30% — 10% of the original one depending on L and (3. 

We observe consistency with chUE except for L = 16, (3 > 6, where the 
physical volume is likely to be too small for chRMT to apply. For the fixed 
point and Neuberger's action an apparent transition to chSE, not observed 
here, took place for small lattices and/or high values of (3 [H]. 

The fit of the chUE distribution provides an estimate for £ which is given 
in Tabled] and in Fig. left (we are using ea = l/y/J3 as abscissa, in the 
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Figure 1: Distribution of the smallest eigenvalue for different lattice sizes 
and P = 4 (quenched case); the continuous curve represents the theoretical 
prediction of chUE. 

ordinate we report = — (ipip) C ont / the equations hold in the weak 

coupling limit of the Schwinger model with one flavor). 

Following RMT the spectrum is to be separated in a fluctuation part 
(conjectured to follow chRMT predictions) and a smooth background. The 
information on the chiral condensate is contained in both. Removing the 
smooth background is done by "unfolding" : the eigenvalues are mapped to a 
variable z, such that the average level spacing is constant. In Fig. |2]we plot 
the microscopic spectral density p s (z) for the unfolded eigenvalues for two 
situations. In the unfolding process the average value of the spectral density 
near A = provides a direct estimate for £ (see (0)). Neglecting finite-size 
effects, this value should agree with the number obtained from the smallest 
eigenvalue distribution in the scaling region of chRMT. We find discrepancies 
between the two determinations for /3 > 4, the "background" value of S (i.e. 
7rp(0)) coming out larger; as a consequence of this discrepancy, the spectral 
density does not follow the predicted shape anymore (Fig. right). We 
interpret this as due to the transition to the full flavor-symmetric situation 
of the continuum limit, which in the quenched case implies trivial doubling 
of the spectrum. 

It has been suggested [TH] that the spectral density of the quenched 
Schwinger model should diverge exponentially with the volume for small 
eigenvalues; our data do not allow us to draw any firm statement concerning 
this. We only observe that for P =2 and 4, L = 24 and 32 give consistent 
values of E, while it is no more the case for the larger value of /3 = 6, this 
indicating a slowing down of the scaling behavior in the weak coupling region. 
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Figure 2: Microscopic spectral density p s (z) for different lattices and val- 
ues of (3 (quenched case); the continuous curves represent the theoretical 
prediction of chUE. 



Topological sector. Fig. El shows the distributions of the first and second 
eigenvalue for configurations with \ v\ = 1. Ideally, the first eigenvalue ought 
to be a zero mode, while the second eigenvalue should follow the predictions 
of chUE for the smallest eigenvalue for \v\ = 1; these may be obtained just 
using the values of £ found in the trivial sector in the formulas for the smallest 
eigenvalue distribution with \u\ = 1 (the continuous line in the Fig. EJ). 

What we see instead is that for small (3 the first eigenvalue follows the 
chUE statistics of the smallest eigenvalue in the trivial sector (dashed curve in 
the figure). Such a behavior was also observed for SU(3), d=4 [21]. Increasing 
/?, the first eigenvalue peaks closer and closer to zero and the second smallest 
eigenvalue approaches the predictions of chRMT for the smallest eigenvalue 
for \v\ — 1 (/3 > 4). We conclude that only towards the continuum limit 
the zero modes correlate clearly with the topological sector of the gauge 
configuration. 

In order to check the systematics of the possible misidentification, we re- 
peated the fits with chRMT shapes without selecting the trivial topological 
sector: even for the smallest value of (3 at our disposal, (3 = 2, we always 
obtain considerably larger values of x 2 /d.o.f. (data with asterisk in Table [TJ 
and the fitting values of £ are inconsistent with the determination in the 
trivial sector. We conclude, that even for moderate values of (3, the iden- 
tification of the topological sector is an essential ingredient for the correct 
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Figure 3: Distribution of the smallest and second (shaded histogram) eigen- 
value for \v\ = 1. The continuous (dashed) curves represent the theoretical 
prediction of chUE for the smallest eigenvalue for \v\ — 1 {y — 0) based on 
the S-value from the fit to the v = data discussed earlier. 
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Figure 4: Results for eigenvalues after re- weighting with the fermion determi- 
nant, corresponding to dynamical fermions. Left: Level spacing distribution 
(unfolded variable), compared with the theoretical expectation for chUEj^j. 
Right: Distribution of the smallest eigenvalue (diamonds with error bars). 
The full and dotted curves are the best fits from chUE with Nf = 1 and 
Nf = 2 respectively, resulting in £ values as given in Table d For compari- 
son we report also the histogram of quenched data (dashed line). 



comparison with chRMT. 

Re-weighting. Of high interest in chRMT is the dynamics of fermions. In 
order to account for dynamical fermions, we re-weight the quenched gauge 
configurations with the fermion determinant; in the case of the statistics 
-P(A m in), for example, this amounts to the following redefinition: 

nAminj ~ EcdeW(C) ' [6) 

where A m i n (C) denotes the smallest eigenvalue for the configuration C. The 
fluctuation of the determinant introduces additional fluctuations in the stud- 
ied quantity (in our case the error typically doubles), which can be traced 
back to the effective loss of statistics of the gauge sample. Loss of statistics 
is also caused by the depletion of the trivial sector (the one under study) on 
large lattices. 
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Figure 5: The condensate £ as obtained from the fits on the smallest eigen- 
value distribution using chRMT predictions, for the quenched (left) and re- 
weighted (right) gauge configurations: L — 16 (boxes), L=24 (diamonds), 
L=32 (triangles); in the dynamical case we plot the fit values obtained from 
comparison with chUE for Nf = 1 (open symbols) and Nf = 2 (full symbols). 

As already discussed in ^3], in the quenched situation we find perfect 
chUE behavior for the level spacing distribution for all parameter values 
studied. This observation persists also for the re-weighted, unquenched data 
(e.g. Fig. HJleft). For the level spacing we always find chUE distribution. 
Such a behavior in a situation of vanishing fermion condensate was found 
already in the Coulomb phase of compact four dimensional QED [22 . 

In Fig. 0] (right) we report the distribution of the smallest eigenvalue for 
L — 16, /3 — 4 and find it clearly different from the quenched distribution. 
As expected from the discussion of Sec. |3] neither the Nf = 1 chUE shape 
nor the Nf = 2 succeed in fitting the distribution. The two different fits 
result in different values of £ reported in Tableland Fig. right. In both 
cases the values of x 2 /d.o.f. of the fit are rather large, confirming that chUE 
does not give the right prediction. The fitted values of £ are well above zero. 
As discussed, strictly speaking £ ought to be zero, a different chRMT and 
universality class applying. In the (true) one flavor situation the continuum 
condensate is — (V'V')cont/ e ~ 0.15989; this situation is described by one copy 
of Overlap or FP fermions [E] . 
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Table 1: Values of the fermion condensate obtained from comparison of lat- 
tice data with chRMT (chUE). In the square brackets we indicate x 2 /d.o.f.. 
Data with asterisk are obtained without selecting a topological sector. 

4 Discussion 

We conclude that for a careful RMT analysis of the eigenvalue spectrum of 
the staggered Dirac operator the identification of the topological sector is 
advisable. Towards the continuum limit the topologically non-trivial sectors 
exhibit an extra peak at small eigenvalues, which we attribute to the "want- 
to-be" zero modes. For a recent discussion on the relevance of quasi-zero 
modes and their chirality in the context of staggered fermions cf. [23 . 

We identified problems due to the general ambiguity of staggered fermions 
between the choice Nf = 1 applying for finite lattice spacing and Nj = 2 d ^ 2 
applying in the continuum limit. For example, the quenched microscopic 
spectral density clearly deviates from chUE for j3 > 4, because the "back- 
ground" estimate of E, 7rp(0), does not match with the "microscopic" one as 
can be obtained e.g. from the smallest eigenvalue distribution. 

Additional unexpected features come into play because of our two-dimen- 
sional context: the Mermin- Wagner- Coleman theorem implies a vanishing 
fermion condensate with the result that chUE fails in the case of dynamical 
fermions. This is indicated by our data, even if the statistics is not good 
enough to draw firm conclusions. This scenario is quite different from that 
found in jH] where the Schwinger model with Ginsparg- Wilson fermions was 
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analyzed; there, genuine single-flavored fermions where described. 
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